Predicted mobility edges in one-dimensional incommensurate optical lattices: An 
exactly solvable model of Anderson localization 
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Localization properties of non-interacting quantum particles in one-dimensional incommensurate 
lattices are investigated with an exponential short-range hopping that is beyond the minimal nearest- 
neighbor tight-binding model. Energy dependent mobility edges are analytically predicted in this 
model and verified with numerical calculations. The results are then mapped to the continuum 
Schrodinger equation, and an approximate analytical expression for the localization phase diagram 
and the energy dependent mobility edges in the ground band is obtained. 
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Anderson localization, the localization of electronic 
Bloch waves due to interference in disordered poten- 
tials, is one of the fundamental quantum phenomena in 
nature and is the transport mechanism behind metal- 
insulator phase transitions in solids pQ. Although this 
mechanism was first proposed over 50 years ago, di- 
rect observation of Anderson localization has been no- 
toriously difficult due to the problems in reliably con- 
trolling disorder in solid-state systems. But recent ad- 
vances in the manipulation of ultra-cold atoms offer a 
completely new, well-controlled tool in directly observ- 
ing such fundamental quantum phenomena. A notable 
example is the recent work done by Billy et. al. who ob- 
served Anderson localization in a diffuse Bose-Einstein 
condensate in a f dimensional (fD) waveguide with a 
disordered laser speckle potential [3J. Another recent 
example is the work done by Roati et. al. who ob- 
served fD Aubry-Andre localization (a phase transition 
closely related to Anderson localization]!}]) of cold-atoms 
in an incommensurate quasi-periodic potential |4j. These 
advances highlight the potential of ultra-cold atoms to 
experimentally probe fundamental quantum localization 
phenomena that previously could only be studied indi- 
rectly or through numerical calculations. Cold atomic 
systems offer precise control of the background potential, 
and the non-interacting limit is easily achievable with a 
very dilute gas of either bosons or fermions. This is the 
context (and the motivation) of the current work, where 
we introduce a new and theoretically exact fD localiza- 
tion model with mobility edges that should be observable 
in cold atomic systems. 

fD localization phenomena are traditionally studied 
with the nearest neighbor tight binding model: 

Eu n = *l(Un-l + Hn+l) + V n U n , (1) 

where t\ is the hopping term representing tunneling be- 
tween nearest neighboring sites and V n is the onsite dis- 
ordered potential [TJ |2] or the incommensurate potential 
[31 H] . The simplicity of allows for exact theoretical 
statements in certain cases. For example, the fD dis- 



ordered Anderson model allows for only localized eigen- 
states at all energies independent of how weak the dis- 
order may be \\. The Aubry-Andre (AA) model with 
the ID incommensurate potential has either all eigen- 
states extended or localized depending on the strength 
of the potential [3] . Thus quantum localization in these 
f D tight-binding models is, in some sense, trivial because 
all states are either localized or extended with no energy 
dependent localization transition as happens, for exam- 
ple, in the 3D Anderson model [5]- However, there is 
growing interest in exploring deviations from the tight- 
binding assumption [SJ |7j . Ultra-cold atoms loaded into 
optical lattices with controllable depths provide an exper- 
imental tool to study transport beyond the tight binding 
regime, where mobility edges are likely. In this Letter, 
we introduce an exact analytically solvable ID localiza- 
tion model which has a energy dependent mobility edge. 
We believe that our model should be realizable in ultra- 
cold ID atomic systems, and we show that our theoretical 
findings extend to the general Schrodinger equation de- 
scription well outside the tight binding regime where cold 
atom localization experiments [5J [3] are typically carried 
out. 

To highlight the new physics that may be observed 
with ultra-cold atoms in shallow lattices, we study local- 
ization in incommensurate lattices with an implicit short 
range rather than nearest neighbor hopping model. In 
particular, we study the tight binding model: 

Eu n = te- pln - n ' l u n ' +Vcos(2tt an + S)u n , (2) 

where a is an irrational number, and p > 0. This is a sim- 
ple exponential hopping generalization of the AA model. 
We prove that this model has energy dependent mobil- 
ity edges (contrary to the AA model), and we verify our 
analytical results by numerically diagonalizing We 
then show that our analytical results from this model can 
be used to predict the energy dependent mobility edges 
in the more fundamental Schrodinger equation for non- 
interacting particles in shallow, incommensurate optical 
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FIG. 1. (Color online) Energy eigenvalues of |2| with 500 
lattice sites and a = (v5 — l)/2 for a) p — 1.5 and b) p = 4. 
The shading of the energy curves indicate the magnitude of 
the inverse participation ratio for the corresponding wave- 
functions. The solid line represents the analytical boundary 
between spatially localized and spatially extended states. 



lattices: 



2m dx 2 



coa(2kr 1 x) + 



Vi 



cos(2kLax))ip{x) — Ei{j(x), 



(3) 



where Vo is the strength of the primary lattice and Vi is 
the strength of the secondary lattice (< Vq). We verify 
these results by numerically solving ^ and examining 
the eigenstates. The Schrodinger equation described by 
Q is often called the two-color or the bichromatic po- 
tential problem 0]. 

The AA model, where the potential in ([!]) is given by 
V n = V cos(2n an + S), has been shown to be self-dual 
under the transformation: 



Ef fD itn{2Tran+8) ifin 



(4) 



a) low E state below edge 






_ b) high E state above edge 




c)state near boundary 




...^ ^ Jli^W) 


ijlkJkL ~ 



50 1DD 150 200 250 3DD 35Q 400 45Q 5DD 

site number 



FIG. 2. Eigenstates of ^ with 500 lattice sites, a = - 
l)/2, V = 1.8, and p = 1.5 for different energy eigenvalues: 
a) low energy localized state below the mobility edge b) high 
energy extended state above the mobility edge c) critical state 
near the mobility edge. 



duality inherent in ([lj does not apply to ^ and in par- 
ticular, the transformation defined by Q does not work 
for the finite range hopping model. 

We now show below that ^ allows for energy depen- 
dent duality points, (i.e a mobility edge). Define the 
parameter p > such that: 



(E + t)- Vcos(27rcm + S) = u 2 T n 



(5) 



cosh(p ) - cos(27rcm + S) 

T„ = ■ s , (6) 

smh(po) 



with tu 2 = ^/(E + t) 2 - V 2 . It follows that [E + t)/V 
cosh(po) and ^ can be cast in the form: 



If we now consider the transformation: 

~ _ im(2Tvan+8)rp 
n 

and note that for p > we have the identity, 



= £ 



p\m \ gim(2TTan-\-5) 



(7) 



(8) 



(9) 



when V = 2ti [3J. Since Q transforms spatially local- 
ized states into spatially extended states and vice-versa, 
it is argued in [3J that all eigenstates are extended for 
V < 2t\ and localized for V > 2t\. The eigenspectrum 
at V = 2t\ has been shown to be singular continuous and 
the eigenvalues produce a Cantor-set structure [8J. This 
sharp transition between localized and extended states 
is often referred to as (AA) duality. Obviously, the AA 



then it follows that the state, u m satisfies the equation: 
u 2 f m u m = J2 te-^ m - m '\u m ,, (10) 



where T m is given by: 

cosh(p) — cos(27ram + S) 



sinh(p) 



(11) 
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FIG. 3. (Color Online) Inverse participation ratios of all 
eigenstates of ([2J with 500 lattice sites and a = (a/5 — l)/2 
for a) p = 1.5 and b) p = 4. The solid curves represent the 
analytical boundary between spatially localized and spatially 
extended states and the dashed lines represent the Aubry- 
Andre condition. 



FIG. 4. Eigenstates of (J3J with 500 lattice sites, a = 
(\/5- l)/2, V = 2E r and Vi = 1.43EV for different en- 
ergy eigenvalues: a) low energy state below mobility edge b) 
high energy state above mobility edge c) critical state near 
the predicted mobility edge. 



We see that ^ is self dual under the transformation ([8]) 
when p — pq. Following AA [3J, we conjecture that all 
states are localized for p > po and extended for p < pq. 
Then it follows that the condition for localization is given 
by the expression: 



cosh(p) 



E + t 
V 



(12) 



It is straightforward to check that this condition becomes 
the AA condition where ^ becomes ([!]) in the limit p — > 
00, (i.e. V = 2ti). The condition given by (12) is the 
central new result of our work, showing that the model 
defined by ^ has an energy dependent mobility edge 
characterized by a transcendental equation. 

To explicitly verify ( 12 ), we numerically diagonalize ([2| 
and study the spatial extent of the wavef unctions. To do 
this we calculate the inverse participation ratio (IPR): 



IPR W = 



(E„K j l 2 ) 



(13) 



where the superscript i denote the i-th eigenstate. The 
IPR of a wavefunction approaches zero for spatially ex- 
tended wavefunctions and is finite for localized wavefunc- 
tions. 

Fig[T]plots energy eigenvalues and the IPR of the corre- 
sponding wavefunctions for ([2| as a function of potential 
strength, V, with a = (a/5 - l)/2 and p = 1.5 (Fig. [TJl) 
and p = 4 (Fig. [TJd) . The solid line in the figure represent 
the boundary given in ( |12[ ). As expected from our con- 
jecture, IPR values are approximately zero for energies 
above the boundary and are finite for energies below the 



boundary, indicating that (12 1 indeed defines the mobil- 
ity edge for pi). 

Fig. [2] displays eigenstates for p — 1.5 and V — 1.8 at 
three different energy eigenvalues which lie above, below 
and near the mobility edge given by (12). We see that 



the wavefunction is localized for low energies (Fig |2ji) , 
extended for high energies (Fig^), and critical near the 
boundary (Fig[2j;). In Fig. [3j we show the IPR of the 
wavefunctions as a function of eigenstate number, i, and 
potential strength, V. The solid curves represents the 



predicted boundary given by (12) and the dashed line 
is the AA self duality condition (V — 2ti). We see in 
Fig. [3] that our mobility edge prediction agrees well with 
the numerical IPR calculations. In the case where p is 
relatively small, the phase diagram shows clear mobility 
edges and differs markedly from the AA condition. For 
large p, the slope of the localization condition is steep 
and is approximately equivalent to the AA condition, as 
expected. 

To understand the relevance of these results to ultra- 
cold atoms in optical lattices, we draw the connection be- 
tween the exponential hopping tight-binding model given 
in ([2| and the fundamental single particle Schrodinger 
equation in ([3]). To do this, we study the ground band 
Wannier functions [5], w n (x) of ^ for Vi = and ap- 
proximate the matrix elements of the Hamiltonian in the 
Wannier basis. Using the Gaussian approximation for 
the ground band Wannier states, the potential strength, 
V, in ([2| is approximated by the expression: 



t 7- V 1 o? 
V ss — exp( , 



(14) 



where E r = (HIcl) /2m is the recoil energy. Also from 
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FIG. 5. (Color Online) Inverse participation ratios of the 
approximate ground band eigenstates of ^ with a = (VE5 — 
l)/2 for a) Vo = 2E r and b) Vo = 5E r . The solid curves 
represent the analytical boundary between spatially localized 
and spatially extended states. 



this approximation, we have for the constant energy dif- 
ference between ^ and (jSj): 



E = (w n \H \w 7 , 



1 



(V e 



+ VVo~E~), (15) 



where Hq is the Hamiltonian corresponding to (|3| with 
Vi = 0. The hopping coefficient, t, can be estimated us- 
ing the deep lattice approximation for the ground band- 
width: 



(16) 



To estimate the exponential decay term, p, we make use 
of Kohn's results on the Kramers' function and its rela- 
tion to the asymptotic behavior of the Wannier functions 
[TUl ITT] . Using the deep lattice approximation for the ef- 
fective mass at the top groundband edge of H$, we obtain 
for p the approximation[T2] : 




p ss cosh 1 ( 1 



2Wn 



), 



(17) 



where Wq is the bandwidth of the ground band and W1/2 
is the width of the first bandgap. The ratio Wi/ 2 /Wq can 
be estimated using properties of the Mathieu functions: 



W 1/2 
W 



-(£)- 
8 E r 



1/4 



exp 




(18) 



Finally, using (17) and (14), the mobility edge for the 
1-D incommensurate lattice Schrodinger equation, ^ , is 
given by: 



2 exp 



{E-E +t) = V 1 1 



2Wn 



(19) 



the Schrodinger equation equivalent of (12) 



To examine the accuracy of (|19j), we numerically in- 
tegrate Q to obtain the energy eigenvalues and wave- 
functions and calculate the IPR (obtained with ( 13 ) by 
replacing the sums with spatial integrals). In our cal- 
culations, we set ki, = 1, a = (\/5 — l)/2, m = 1. 
The size of the system is given by L = Na where a 
is the lattice constant. N is chosen to be 500 and ^ 
is sampled over 80,000 points. Fig. [3] gives eigenstates 
at three different energy eigenvalues for Vo = 2E r and 
V± = lA3E r (similar to Fig [2]). Similar to the results in 
the tight binding model, we see that for a fixed potential 
strength an eigenstate can be localized for low energies 
(Fig[4ji), extended for high energies (Fig [4^), and critical 
near the boundary (Fig [4}:). Fig. [5] gives calculated IPR 
values as a function of eigenstate number and V\ for the 
first N eigenstates (equivalent to the ground band when 
Vi = 0) The solid curves give the analytical boundary 
between localized and extended states as given by ( 19 1 



We see in Fig [5] that our analytical prediction is in good 
agreement with our IPR calculations. We also note that 



( 19 1 is dependent on incommensuration and may predict 
no localization transition for a 2 / yJV^/ E r ^> 1, where the 
slope of the boundary in E — Vi space is essentially flat. 
This is consistent with numerical results reported earlier 
in [6] where localization transitions in |3]) are observed to 
be dependent on incommensuration for shallow lattices. 

We have predicted the existence of an exact, analytic, 
energy dependent mobility edge in an incommensurate 
ID model, which should be experimentally accessible in 
cold atomic systems. This mobility edge is the energy de- 
pendent generalization of the AA self-duality concept and 
as such, the eigenstates precisely at the mobility edge are 
critical, (i.e. neither localized nor extended) with the mo- 
bility edge spectrum being singular continuous. The fact 
that our predicted tight-binding mobility edge survives 
the continuum Schrodinger equation limit indicates that 
our exact prediction should be observable in cold atomic 
systems. We note that we have explicitly numerically 
verified, |13j . that the exponential hopping constraint is 
not crucial - in fact, any bichromatic Schrodinger equa- 
tion with hopping amplitudes falling off in some manner 
(e.g. power law, Gaussian, etc.) would approximately 
exhibit the predicted ID mobility edges. We also note 
that the slowly varying trap potential confining the cold 
atoms, not included in our calculations, should not affect 
our conclusions because it only acts as a finite size cut 
off for the extended states, which is incorporated in our 
numerical simulations. 
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